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ABSTRACT 

The frequencies of kHz QPOs are widely interpreted as being indicative of the values of charac- 
teristic frequencies related to orbital motion around neutron stars, e.g., the radial epicyclic frequency. 
In regions directly exposed to the radiation from the luminous neutron star these frequencies change 
with the luminosity. Including radiative corrections will change the neutron star mass value inferred 
from the QPO frequencies. Radiative forces may also be behind the puzzling phenomenon of parallel 
tracks. 
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1. Introduction 

"Twin peak" QPOs with a pair of physically connected frequencies Vu , Vl 
are observed in several X-ray binaries (e.g., van der Klis 2000). Physical mod- 
els involve characteristic frequencies related to orbital motion, either directly (e.g., 
Morsink & Stella 1999), or indirectly through their effect on disk oscillation prop- 
erties (e.g., Wagoner 1999, Kato 2001). The values of these frequencies depend on 
the mass of the central neutron star or black hole. 

The frequency ratio in all black hole sources in which the twin peak QPOs have 
been detected is Vu/vl = 3/2 exactly. In several neutron star X-ray binaries the 
ratio has approximately the same value, Vu/vl ~ 3/2. Kluzniak & Abramowicz 
(2000), Abramowicz & Kluzniak (2001), and Kluzniak & Abramowicz (2002), 
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proposed that the twin peak QPOs may be explained by a non-linear 3:2 resonance 
in the epicyclic oscillations of accretion disks around the compact object in these 
sources. The resonance occurs because of properties of strong gravity in general 
relativity and hence the model predicts that the QPO frequencies scale inversely 
with the gravitational mass of the object, 

v~l. (1) 
M 

However, Abramowicz, Bulik, Bursa & Kluzniak (2003) pointed out that when 
this model is applied to the atoll and Z sources, where Vu ~ 900 Hz, a rather low 
value for the neutron star mass is obtained M < 1.2Mq. In the case of Sco X- 
1 those authors found that the epicyclic resonance model predicts the mass to be 
M = I Mq . A similar problem occurs for most of the neutron star sources, as 
exhaustively discussed in a recent paper by Urbanec et al. (2010). Note that the 
black hole QPOs have lower frequencies, and accordingly the implied masses can 
be much higher, and for particular values of the Kerr parameter can be made to 
correspond to the measured black hole masses. 

Discussing this problem, Urbanec et al. (2010) concluded that the epicyclic 
resonance model cannot explain the QPO phenomenon in neutron star sources, 
unless the oscillating accretion disks around neutron stars experience some "fairly 
non-geodesic" effects, i.e., when there are strong non-gravitational forces acting on 
the disk matter. 

In this Note we suggest a natural origin for such "fairly non-geodesic" effects. 
We expand on an idea already suggested by Yu (2008, 2011) that radiation force 
affects the orbital frequencies in QPO sources, and show (in Section 3) that when 
the illumination of the accretion disk by the radiation emitted by the neutron star 
is properly included, the 3:2 epicyclic resonance model implies higher neutron star 
masses. The main point is that radiation forces reduce the effective gravity. The 
observed frequencies Vu ~ 1 kHz may then be consistent with the canonical mass 
value ~ IAMq , and even higher ones. We note that Orosz & Kuulkers (1999) 
measured a mass of 1.78 ±O.23M0 for the neutron star Cyg X-2, which is known 
to be a high luminosity twin kHz QPO source (Wijnands et al. 1998). 

2. Illuminated epicyclic oscillator 

Figure [Uillustrates the effect discussed here, for the epicyclic resonance model. 
The interior of a geometrically thin accretion disks around a neutron star is affected 
by radiation of the luminous boundary layer only in a very small region near the 
inner edge rjn. This is because the optical depth in the radial direction is very 
large, and the disk is self-shadowing. This has two obvious consequences. Firstly, 
the radiative drag (Poynting-Robertson effect) discussed by Walker & Meszaros 
(1989), and others, cannot affect motion of matter inside the disk at the location 
of the 3:2 resonance radius r3-2 ~ 2rin. Secondly, the vertically oscillating part 
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Fig. 1. Oscillating accretion disk illuminated by the neutron star boundary layer. The inner disk is 
optically thin in the vertical direction, and optically thick in the radial direction (note the shadow in 
the cut-away). In the figure, the disk has been truncated for clarity. 



of the disk (shown in Figure [T] as a toroidal ring) periodically emerges from the 
disk shadow, and hence, is illuminated by the boundary layer. If the oscillating 
torus is optically thin (as will be the case when the disk is optically thin in the 
vertical direction) the emergent ring will briefly be subjected to a radiative force 
proportional to the mass of the ring that is momentarily outside the shadow of the 
disk. Thus, the effective gravity on the ring is reduced by a factor proportional to 
the luminosity of the boundary layer. 

In passing, we note that the radiation force caused by this illumination is a func- 
tion of the vertical displacement, 7=7 (8z) . This leads to a coupling between the 
vertical and radial oscillations that can be described by a forced oscillator equation, 

5r + (o25r=ir,(|5z|), (2) 

where CO,- is the radial epicyclic frequency corresponding to pure gravity (absence 
of illumination), and the form of the forcing function is to be calculated from the 
exact geometry of illumination. Such a coupling is a desirable feature in the QPO 
resonance model. 
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3. Corrections to the mass estimated from tlie oscillatory frequencies 

Observations of the response of the twin kHz QPOs to the flux of the normal 
branch oscillations (NBOs) in the manner of an anti-correlation between the upper 
kHz QPO frequency and the NBO flux in Sco X-1 provided evidence that radiation 
force affects the frequency of the orbital signal on sub-second time scales (Yu, van 
der Klis & Jonker 2001). More solid evidence is found for this anti-correlation in 
the atoll source 4U 1608-52 (Yu & van der Khs 2002). Yu (2008 and 2011) sug- 
gested that the radiation force can cause the apparent QPO coherence difference be- 
tween Z and atoll sources and the apparent coherence drop in several well-observed 
atoll sources (Barret et al. 2005, 2008), and that the radiation force effect can be 
used to constrain neutron star masses and radii if both radiation force and orbital 
frequency can be measured. For a clear positive correlation between luminosity 
and QPO frequency see e.g.. Fig. 1 of van der Klis (2001). 

The radiation force effects that are suggested by the correlated behavior of the 
kHz QPO frequency and the X-ray flux in both atoll sources and Z sources provide a 
natural origin of the "fairly non-geodesic" effects invoked by Urbanec et al. (2010). 
This is because radiation leads to a reduced "effective gravity" and consequently 
the mass M that appears in ([T]) is not identical with the true gravitational mass of 
the star M^, , but instead it is the smaller "effective mass" 

M = M, fl-£-^y (3) 

Here, L* is that part of the total observed luminosity L that influences the disk 
oscillations, and the part that does not will be Li = L — L^, while L^dd is the Ed- 
dington luminosity, and £ < 1 is a factor that depends on details of the radiation 
field geometry and of the radiation transfer in the accretion disk. Accordingly, the 
mass of the neutron star predicted by the epicyclic model when radiative corrections 
are included is. 



Me 



fy.2 



^Edd 



-1 



(4) 



Here fj,-2 ~ 900 Hz for the 3:2 resonance in Schwarzschild metric. For different 
orbital resonances (discussed by several authors) the value of this coefficient will 
be different. 

One may estimate values of £ and L, from observations. From @ it follows, 
in the case of small zL^, , that 

/ £Mo \ Me f eLi \ 



If the values of £ and L, do not change much during the observation period, 
and L, changes from one observation to another (indicated by the / label), then 
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V = aL + bi. In this case one may estimate £ and L, from the observable parame- 
ters a = a{e,M^), and bi = bi{e,M^,Li) . Thus, Eq. ^ may explain the observed 
parallel tracks (Mendez et al. 1999), in which the QPO frequency is proportional 
to the luminosity on short timescales, but for successive observations the lines on a 
luminosity-frequency plot are displaced in frequency by different offsets. 

4. Conclusions 

In conclusion we want to stress three points: 

• The radiation force from the luminous neutron star will affect the orbital, and 
related frequencies, most strongly at the inner disk edge. 

• In all previous estimates of the neutron star mass based on the 3:2 epicyclic 
resonance QPO model it was assumed that £ = 0, and this implied masses 
that were too low. The inclusion of radiation force effects will lead in a nat- 
ural manner to an increased mass prediction. The brighter the neutron star 
source is, the more the inferred mass will be raised, since the mass was un- 
derestimated in the original models by an amount proportional to the strength 
of the radiation force effect. For Sco X-1, if eL = O-SLndd. the mass estimate 
would double, yielding a mass up to 2 solar masses. For atoll sources, the 
mass estimate would be raised by 10-20%, e.g., by 0.1 — O.4M0, assuming 
L is about 10-20% Lndd- 

• The inclusion of the radiation force may provide an explanation for the par- 
allel tracks phenomenon. 
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